Large scale simulations of the jet-IGM interaction 
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Abstract. In a parameter study extending to jet densities of 10 _ times the ambient 
one, I have recently shown that light large scale jets start their lives in a spherical bow 
"^J ' shock phase. This allows an easy description of the sideways bow shock propagation 

in that phase. Here, I present new, bipolar, simulations of very light jets in 2.5D 
and 3D, reaching the observationally relevant scale of > 200 jet radii. Deviations 
from the early bow shock propagation law are expected because of various effects. 
The net effect is, however, shown to remain small. I calculate the X-ray appearance 
of the shocked cluster gas and compare it to Cygnus A and 3C 317. Rings, bright 
spots and enhancements inside the radio cocoon may be explained. 
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1. Introduction 

o . 

^£ ■ X-ray studies of galaxy cluster centers containing a radio jet have 

shown that the jets have a considerable impact on the cluster gas 
(Carilli et al., 1994; Carilli et al., 2002; Blanton et al., 2001; Smith 
et al., 2002), forming rings, apparent spirals, and aligned features. 
Such systems have been claimed to be associated with very light jets 
(Clarke et al., 1997; Rosen et al., 1999; Krause and Camenzind, 2003). 
Ci ■ Parameter studies of very light jets have been carried out recently (Car- 

valho and O'Dea, 2002a; Carvalho and O'Dea, 2002b; Saxton et al., 
2002; Krause, 2003; Zanni et al., 2003). There, it has been shown that 
very light jets first form spherically symmetric bow shocks (Krause, 
' 2003). In that phase, they follow the expansion law (derived for a strong 

bow shock, which is applicable here): 

M(r')r' dr' = 2 /* dt' /* E(t")dt" , (1) 
o Jo Jo 

where Jv[(r) is an arbitrary spherically symmetric ambient gas mass 
distribution and E(t) is the energy injection law. Here, I show 3D and 
2.5D large scale simulations, updating the X-ray appearance of the 
shocked IGM, and exploring the accuracy of the spherical expansion 
law for the bow shock at late times. 
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Number density at t= 2.04 Myr 
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Figure 1. Logarithm of number density (slice) for the 3D run at 2.04 Myr. 



2. 3D Simulation 



2.1. Simulation Setup 



A cylindrical grid was used for the jet simulation (compare Krause 
and Camenzind (2002)). The size of the computational domain was: 
Z £ [-69 kpc,69 kpc], R £ [0,57 kpc] and <j> £ [0,2vr]. 2042, 805, 
and 57 grid points were used in the Z,R and <f> directions, respectively. 
The full extend was not reached because of long computation times 
and breakdown of the supercomputer. With a jet radius of rj = 0.55 
kpc, this gives a resolution of 8 points per beam radius (ppb). The grid 
was initialised with an isothermal King cluster atmosphere: p c {r) = 

( 2 \ -3/3/2 

p e ,o (1 + ^z) j where r = yi? 2 + Z 2 denotes the spherical radius, 

p Ci o = 1-2 x 10 -25 g/cm 3 is the characteristic density, j3 = 0.75 and 
a = 35 kpc is the core radius. In order to break the bipolar and 
axial symmetry, this density profile was modified by random pertur- 
bations. The jet is injected in the middle of the grid in the region 
Z G [—0.55, 0.55 kpc] ,#£[0,0.55 kpc] , and </> e [0, 2vr] . This region has 
the constant values: p ]e t = 6.68 x 10 -28 g/cm 3 , vz = ±0.4c, c being the 
speed of light. The kinetic jet luminosity is L^ n = 1.04 x 10 46 erg/s for 
both jets together. The pressure was set in order to match the external 
pressure at that position. This gives a slightly varying density contrast 
across the grid of r\ = p ]e t/Pext ~ 7 x 10 -3 and an internal Mach 
number M = 10. The temperature in the external medium is set to 
3 x 10 7 K. The cooling time in the shocked cluster gas is approximately 
100 Myr. The jet is expected to propagate through the whole volume 
in 10 Myr. So, cooling by bremsstrahlung marginally influences the 
state of the gas. This was taken into account. In order to keep the 
system in hydrostatic equilibrium, gravity by an assumed dark matter 
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Bow shock shapes at t = 2.04 Myr: elliptical + cigar 




Z [kpc] 

Figure 2. Bow shock shapes for the 3D run, compared to ellipse and parabola. 

distribution had to be applied. The non-relativistic code NIRVANA (?) 
was used for the computation. 

2.2. Results 

The final snapshot at t = 2.04 Myr is shown in Fig. 1. The cocoon is 
nicely placed around the jet beam. Shear instabilities show up promi- 
nently. They still grow towards the center, and develop into long fingers 
at the innermost positions. The pressure shows a regular spacing of 
shock compressed and rarefaction zones in the beam. High pressure 
regions are small and show up only at the end of the beams, where the 
Mach disk is located. The central region, with a diameter of roughly 
10 kpc, is now dynamically calm. No large Mach numbers are observed 
there, and the pressure is approximately constant. 

2.2.1. The shape of the bow shock 

Figure 2 shows the bow shock shape for the final snapshot in detail. 
It has an axisymmetric part in the middle, where it can be well repre- 
sented by an ellipse. The elliptical shape ends at \Z\ ~ 10 where two 
cigar like extensions join the bow shock. They are 3D in nature, and 
can be represented, on average, by a parabola of rank three. 

2.2.2. Emission maps 

The emission due to bremsstrahlung was integrated for different view- 
ing angles (see Fig. 3). The general X-ray emission properties of shocked 
ambient gas have been discussed by Clarke et al. (1997), which has 
been updated recently by Zanni et al. (2003). The idea is that the 
gas is pushed aside by the jet cocoon. Depending on its compression, it 
may form X-ray deficits at the location of the cocoon, and bright shells 
at the edges. The critical parameter is the relative shell thickness, £, 
defined as the width of the shocked ambient gas region divided by the 
local bow shock radius. X-ray deficits could be observed for £ > 38%, 
for sources at inclination i = 90°. Here, £ is comparatively low. Hence, 
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Figure 3. Brcmsstrahhmg emission maps for the 3D run at t = 2.04 Myr. 
From top to bottom, the viewing angle is 0°, 10°, 30°, 60°, and 90°. The left 
column shows the emission map, the middle one a vertical, and the right one 
a horizontal slice through the center. The undisturbed emission is indicated. 



at high i the X-ray surface brightness never falls below that of the 
undisturbed King atmosphere, but the deficit is pronounced for low i 
(compare Fig. 3). 

The two phases of the bow shock, cigar and elliptical (see sect. 2.2.1), 
show up prominently in the emission maps. They form circular and el- 
liptical rings, depending on the viewing angle. Where the rings partially 
overlap, they are brighter, producing the impression of ring segments 
(e.g. Fig. 3, 10°). The structures can also intersect on the line of sight, 
producing bright spots (Fig. 3, 30°). The pole on figures show at least 
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two rings: one from the cigar phase and one from the inner elliptical 
part. 



3. 2.5D simulation 

3.1. Simulation Setup 

In order to study the jet evolution on larger scale, an axysymmetric 
(2.5D) simulation was performed with initial conditions similar to the 
3D simulation in the previous section. The 20 ppb simulation was run 
for 20 Myrs of simulation time. During that time the jet reached an 
extent of 110 kpc which corresponds to 220 jet radii. The King atmo- 
sphere parameters are: p e> o = m p /cm -3 , a = 10 kpc, (3 = 0.75, and 
T = 3 x 10 7 K. The jet is injected with a density of pj c t = 10 -4 x p ej o, 
the sound speed in the jet was set to 0.2c, and the jet's Mach number 
to M = 3. 

3.2. Results 

Logarithmic density and integrated X-ray emission after 20 Myr are 
presented in Fig. 4. Krause (2003) could only reach an evolution to up 
to 5 Myr. In this early phase, the bow shock is spherical, its radius 
following an expansion law given by (1). At later times the cocoon 
transforms via a conical phase towards a cylindrical one (Fig. 4). 

3.2.1. Pressure evolution 

The average jet pressure is the driving force of the inner elliptically 
shaped part of the bow shock. The pressure in the jet system mono- 
tonically decreases with radius. Close to the axis, the pressure is higher 
because of shocks in the beam region. In the shocked ambient gas 
region, a new equilibrium of gravity and pressure appears. The smallest 
pressure values are located at the bow shock, roughly 20% below the 
average. 

It has been found that the sideways expansion of the bow shock fol- 
lows the blastwave's equation of motion (Krause, 2003), in the spherical 
phase. The accuracy of this law will be checked in the following also for 
the larger 2.5D simulation. In this case, the bow shock has propagated 
more than three core radii in the sideways direction. In the spherical 
approximation, the average pressure inside of the bow shock is given 
by (neglecting energy stored in the beam): 

Pj = ( 7 -l)(Lt-^ 2 /2)/^, (2) 
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X— ray emission at t=20 Myr 




Number density at t= 20 Myr 
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Figure 4- Final snapshots of the 2.5D simulation. Top: Line-of-sight integrated X-ray 
emission due to bremsstrahlung. Bottom: Logarithm of the number density. 



where Vj is the jet volume (everything inside the bow shock). The power 
L includes all sources of energy, i.e. the flux of kinetic and internal 
energy through the jet channel, the flux of internal energy entering 
through the surface of the bow shock, and the energy lost by work 
against the gravitational field. Then, (2) can be evaluated, where v 
and t are given by (1): 

t j.2 / o j-r \ 1/3 

v= t= v / M(r')r'dr') . (3) 

Mr \L Jo J 

Figure 5 shows this analytical estimate together with the data from 
the simulation. Here, r was related to time via measurement from the 
simulation. The agreement is quite good, in general. The analytical 
formula follows the slope of the simulation data, but underestimates it 
by up to » 20%. 
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Table I. Fit parameters for the bow shock position in the 2.5D simulation. 
The star denotes a fit with fixed a — 0. Two stars denote fits to the spherical 
approximation with fixed a — 0. 



time range [Myr] 


a 
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time range [Myr] 


a 


b 


c 


[0:5] 


2.00 


2.67 


0.76 










[0 : 5f 





4.66 


0.35 


[15 : 20]* 





2.99 


0.78 


[0 : 5]" 





3.22 


0.67 


[15 : 20]" 





2.36 


0.81 



3.2.2. Sideways motion of the inner bow shock part 
The bow shock propagation was locally fitted by a function of type 
a + bt c . The resulting parameters for the different regions are shown 
in Table I. Usually, a = (fixed), since only the late evolution of the 
jet is studied. Only for the time span up to 5 Myr a fit with o^0 has 
been included because here it is possible that effects from the initial 
conditions still dominate the propagation. For comparison, also fits to 
the detailed spherical approximation are given, computed by applica- 
tion of (3). According to that, the exponent for the first five million 
years should be 0.67. Using the pure power law, an exponent of 0.35 is 
achieved in the simulation data. Allowing for the radial offset gives a 
best fit exponent of 0.76. Since the exponent of 0.35 is much below any 
expectation, it follows that the initial conditions are still important in 
that phase, and the fit with offset is more appropriate. The concurrence 
of the curves increases with time. For the last five million years, the 
exponent for the spherical approximation (0.81) exceeds the one for the 
power law fit of the simulation data (0.78). From the increasing aspect 
ratio, an exponent lower than the one of the spherical approximation 
was expected. The simulation shows that the effect is small. 

4. Discussion 

The bipolar simulations in connection with detailed studies of the early 
lives of very light jets reveal two parts of the bow shock: an inner 
elliptical part, and an outer cigar shaped one. These parts also appear 
in the X-ray emission. When viewed from an appropriate angle, they 
can appear as partial rings and bright spots. I suggest that the two 
X-ray rings in 3C 317 (Blanton et al., 2001) are caused by this effect. 

The X-ray appearance of the 2.5D simulation reproduces some im- 
portant details of Cygnus A's X-ray emission (Smith et al., 2002): 
the elliptical deformation in the shocked ambient gas region, the bright 
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Figure 5. Left: Average jet pressure over simulation time. The stars show the values 
measured in the simulation, crosses mark the pressure according to a spherical 
approximation, plusses show the relative difference of the former in percent. Corre- 
sponding symbols are connected with solid lines. The lines show fits to the pressure 
measured from the simulation. The fits are: 32.84 1~ 69 (0-5 Myr, solid line), and 
227.95 i -1,66 (15-20 Myr, dashed line). The best fit for the spherical approximation 
in the range 15-20 Myr is: 126.05 1 -1 ' 50 (not shown). Right: Bow shock radius at Z=0 
versus time (squares) with fits and compared to spherical approximation, including 
all power sources (see text, plus-signs). The three fits are: 4.66t 0,53 (0-5 Myr), 
2.99t ' 78 (15-20 Myr), 2.36 1 81 (15-20 Myr, spherical approximation). 



X-ray filaments inside the cocoon, and the fork like structure around 
the cocoon. The bow shock can be located at the interface, where the 
elliptical X-ray contours meet the spherical ones from the unaffected 
cluster gas. This good agreement also supports the idea of a very light, 
relativistic, and magnetised jet in Cygnus A. 
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